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Abstract 

We revisit calculations of the cosmogenic production rates for several long-lived iso- 
topes that are potential sources of background in searching for rare physics processes 
such as the detection of dark matter and neutrinoless double-beta decay. Using up- 
dated cosmic-ray neutron flux measurements, we use TALYS 1.0 to investigate the 
cosmogenic activation of stable isotopes of several detector targets and find that the 
cosmogenic isotopes produced inside the target materials and cryostat can result in 
large backgrounds for dark matter searches and neutrinoless double-beta decay. We 
use previously published low-background HPGe data to constrain the production of 
"^H on the surface and the upper limit is consistent with our calculation. We note 
that cosmogenic production of several isotopes in various targets can generate po- 
tential backgrounds for dark matter detection and neutrinoless double-beta decay 
with a massive detector, thus great care should be taken to limit and/or deal with 
the cosmogenic activation of the targets. 
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1 Introduction 



The direct detection of weakly interacting massive particles (WIMP) [T] and 
the search for neutrinoless double-beta (Oz//?/3) decay [2] are fundamentally im- 
portant for physics beyond the standard model of particle physics. The direct 
detection of WlMPs would help to determine the mass and cross section of a 
WIMP, while any observation of neutrinoless double-beta decay would clearly 
show that the neutrino is a Majorana particle and that the lepton number 
is not conserved. Both types of experiments are searches for extremely rare 
signals and thus require large mass exposures of sensitive detectors with small 
internal backgrounds and sufficient shielding against external backgrounds at 
a deep underground site. 

The physics goals of upcoming double-beta decay experiments are to probe 
the quasi-degenerate neutrino mass region as low as 100 meV and demonstrate 
that backgrounds can be achieved at or below 1 count /ton/year in the Oiy(3P 
decay peak region of interest (ROI). To realize these goals, these experiments 
must construct detectors in an ultra- low background structure [3]. Although 
these experiments are foremost neutrino mass experiments, they may poten- 
tially contribute to dark matter searches [IfS] . In contrast, dark matter exper- 
iments will not necessarily be built of a material composed of Ou^P isotope. 
In any case, cosmic-ray produced long-lived isotopes are potential sources of 
background for either type of experiment. 

In this paper, we revisit old calculations of the production rates for critical 
cosmic-ray produced isotopes. The need for this revisit arises, because of recent 
improvements in the understanding of the cosmic-ray neutron flux and, in the 
case of '^H, the availability of reaction codes that fully identify all the reaction 
products in the final state. In Section [2] we evaluate, several long-lived isotopes 
that are produced in germanium by cosmic ray neutrons while the material 
resides on the Earth's surface. We then validate our calculation by comparing 
them to the measured production rates in Section [31 We discuss the effect 
of cosmogenic activation of natural xenon and various other targets on dark 
matter detection in Section HI Finally, we give our conclusions in Section [5l 



2 Cosmogenic production in germanium 

2. 1 Cosmogenic production of radioactive isotopes 

The early work of Avignone et al. [13] showed that isotopes produced in Ge 
by fast cosmic-ray neutrons could create background in double beta decay 



2 



and dark matter experiments. Both the Heidelberg-Moscow [6] and IGEX [7] 
experiments observed 7-rays from such isotopes (e.g. ^^Ge, ^"Co). Thus, mo- 
tivated the need to hmit the cosmogenic activation, a better understanding of 
the production rate of cosmogenic isotopes is required. The production rate, 
Ri, of the radioactive isotope i can be calculated according to 



where Nj is the number of target nuclear isotope j, aij is the neutron excitation 
function for the product i from target j, and (p is the cosmic neutron flux. We 
ignore the contribution of cosmic protons because the flux is much smaller. 
Ziegler carried out a comprehensive study on neutron cosmic ray flux [S] and 
pointed out that some of the data from early measurements PITU] is incorrect 
or of marginal quality [11]. Several cosmic neutron fluxes [9IT2] were used 
in the evaluation of cosmogenic production in both natural germanium and 
enriched ^^Ge conducted by Avignone et al. [13] while significant differences 
emerged in the production rate of ^^Ge and ^''Co due to the varying values for 
the cosmic neutron flux. The quality of the neutron flux data has significantly 
improved due to new measurements [H]. Improved recent measurements by 
Gordon et al. [15] show that the flux density spectrum at sea level can be 
parametrized as 

^{E) = 1.006 X 10-6e-0-351n^i5+2.14511ni? ^ ^_q^^ ^ ^Q~3g-0.4106 In^ i?-0.6671n (3) 

where E is neutron kinetic energy in MeV and (j) in units of cm~^s~^MeV~^. 
This parametrization function agrees with the data within ~2% accuracy as 
shown in Fig. [H 

Note that the variation of the cosmic ray muon and neutron fluxes at vari- 
ous locations around world have been reported by Ziegler [S] and Gordon et 
al. [15]. Cosmic ray muon and neutron fluxes depend strongly on the altitude 
and the variation is described as a function of altitude in Ziegler's paper [8]. 
In addition, the variation from different locations at the sea level caused by 
geomagnetic rigidity is also substantial [51IT5] . In the northern hemisphere, 
this variation is within 10% [12]. This work applies the measured cosmic ray 
neutron flux in the northern hemisphere described in Eq.([2]). 

We use the TALYS code 1.0 [|16] to generate the excitation functions of isotopes 
produced by neutrons on stable isotopes of germanium and copper. TALYS is 
a nuclear reaction program which simulates nuclear reactions that involve neu- 
trons, photons, protons, deuterons, tritons, ^He- and alpha-particles, for target 
nuclides of mass 12 and heavier. In the TALYS code 1.0, a suite of nuclear- 
reaction models has been implemented into a single code system which enables 
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Fig. 2. Co production cross section as a function of neutron kinetic energy. 

to evaluate nuclear reactions from the unresolved-resonance region up to inter- 
mediate energies. We show, for example, the excitation functions of ^^Co and 
by neutrons on germanium isotopes in Fig. [2] and [31 respectively. Based on 
Eq. ([T]), the production rates of cosmogenic isotopes in natural and enriched 
(nominally 86% ^^Ge and 14% '^^Ge) germanium are estimated and tabulated 
in Table [TJ Isotopes with half-lives above 10,000 years are not listed because 
their contribution to the backgrounds are expected to be negligible due to 
their low decay rate. The production rates of ^°Co and ^^Ge in germanium 
agree within a factor of 2 with Ref. [TTfTS] , while the rate of is a factor 
of 7 smaller than the estimates in Ref. [T31IT9] . where the old cosmic neutron 
flux from Ref. [91IT2] or [20] is used. Compared to the old neutron-flux mea- 
surements, the recently measured flux is smaller at energies below 50 MeV, 
but larger between 50 MeV and 1 GeV. ^if can produce background for dark 
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Fig. 3. H production cross section as a function of neutron kinetic energy, 
matter detection via /3-decay with an end point at 18.6 keV. 
Table 1 

The calculated production rates (per day per kg) of cosmogenic isotopes in natural 
and enriched (86% ^^Ge and 14% ^^Ge) germanium using TAYLS and the neutron 
flux model described in the text. Very long and short lived isotopes are not listed 
because their contribution to the backgrounds are expected to be negligible. Also 
shown are isotopes which are produced by cosmogenic activation of copper. 



Cosmogenic Isotope 


Production Rate (/ (k 


Z day)) 




Natural Ge 


Enriched Ge 


Natural Gu 


h/2 




41.3 


7.2 




270.8 d 




2.0 


1.6 


46.4 


5.2714 y 




13.5 


6.7 


56.2 


271.79 d 


55 


8.6 


3.4 


30.7 


2.73 y 




2.7 


0.87 


16.2 


312.3 d 




37.1 


20.0 




244.26 d 




1.9 


1.8 




100.1 y 




27.7 


24.0 




12.33 y 



2.2 An Upper Limit on the Tritium Production Rate in Enriched Ge 



For the low-energy data analysis of Ge-detector double-beta decay experi- 
ments, the question of the tritium content is an important issue. The Heidelberg- 



Moscow [21] and IGEX [22] experiments have both pubhshed low-background, 
low-energy spectra from Ge detectors operated deep underground. The Heidelberg- 
Moscow data has a low-energy threshold of 9 keV, whereas the IGEX data 
reaches lower to 4 keV. Figure H] shows the IGEX spectrum. Overlaid on that 
spectrum is the tritium /?-decay spectrum normalized to 250 counts with a 
constant background level of 2.5 counts/keV added. The data is 80 kg-d of 
exposure. The " fit" was done by eye: that is the normalization of the tritium 
curve was determined by adjusting it until it passed approximately through 
the measured data points. 




10 20 30 

Energy (keV) 



Fig. 4. The low-energy IGEX spectrum (80 kg-d) |22j with the fit-by-eye to a tritium 
spectrum. 

If one assumes that the normalized curve accurately represents the data, then 
the assigned 250 counts represents the maximum number of tritium events 
that can be contained within the data. Since it is likely that other sources of a 
signal (noise, low energy x-rays, etc.) are present, we consider the 250 counts 
to be an upper limit on the spectral contribution due to tritium. 

Tritium will be produced within the Ge by high-energy neutron-induced in- 
teractions while it is stored above ground. The Ge detector was underground 
for about 1 year prior to this data being obtained and the enriched Ge was 
above ground 3-5 years after enrichment and prior to going underground |23j . 
Therefore the tritium was produced for 3-5 years and then decayed away for 
a year before counting. 

To convert the 250 counts to an upper limit on the tritium production we use 
3 years for the exposure time as this gives the largest production rate. The 
tritium production rate (k) depends on the number of counts (dN) observed 
during the counting time (dt), the exposure time to cosmic rays (texp) and the 
cool-off time before counting starts (tdec)- The production rate is then given 
by: 
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j^_dN 1 

"^(l-exp^)exp^ 
^ ^ 250 1 

<80kg-d(l-exp^)exp^) 
k <21 tritium atoms/kg-d 

where the mean hfe of tritium is 17.79 y. This resulting upper limit on the tri- 
tium production rate in Ge due to cosmic rays (21 tritium atoms/kg-d) is con- 
sistent with our TALYS code calculation of 24 tritium atoms/kg-d. The rate 
is much lower than a previous calculation [13] of 110-140 tritium atoms/kg- 
d, which couldn't discern the fraction of events of ^^Ge(n,x)^°Ge that led to 
tritium production. 



3 Comparison between this work and the previous measured rates 

Avignone et al. have measured and calculated cosmogenic production rates 
for some isotopes utilizing natural germanium detectors [13]. Their results are 
compared to this work in Table O A reasonable agreement can be found for 
several isotopes except for . The differences in the rates of production 
is due to the difference in the cosmic neutron energy spectra applied in two 
calculations. 

Table 2 



The production rate of isotopes in natural germanium. 



Cosmogenic 
Isotope 


Measured Rate 
(/(kg day)) m 


Calculated Rate 
(Hess Model [9]) 
(/(kg day)) m 


Calculated Rate 
(Lai Model [12]) 
(/(kg day)) m 


Calculate Rate 
(This work) 
(/(kg day)) 






~210 


~178 


27.7 




3.3±0.8 


2.7 


0.93 


2.7 




38.0±6.0 


34.4 


24.6 


37.1 


68Ge 


30±7 


29.6 


22.9 


41.3 



It is worthwhile to mention that our results agree within a factor of two with 
a recent calculation by Barabanov et al. [23]. Table [3] shows the comparison. 
The difference between two calculations are mainly due to the use of different 
cosmic ray flux values. 



(3) 

(4) 
(5) 
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Table 3 

The production rate of isotopes in germanium 



Cosmogenic Isotope 


Natural Ge (/(kg day)) 


Enriched Ge (/(kg day)) 




Ref. [M] 


This work 


Ref. [21] 


This work 




2.86 


2.0 


3.31 


1.6 


68Ge 


82.7 


41.3 


4.32 


7.2 



4 Cosmogenic production in natural Xe and other targets 



Liquid noble gases, such as liquid xenon (LXe) [25|26II27] and liquid argon 
(LAr) [28|29|30j . have shown excellent pulse shape discrimination capabili- 
ties. Liquid cryogens offer the possibility to construct ton-scale target mass 
detectors PT]|32|33] at a reasonable cost. We list in Table H] the cosmogenic 
production rates of isotopes which can produce potential backgrounds for 
xenon-based dark matter detection experiments. Because the interaction rate 
of WIMPs in xenon decreases dramatically with the nuclear recoil energy, the 
sensitive energy region for xenon-based dark matter experiments is the very 
low energy region. Thus, /9-decay with the end-point energy of 18.6 keV 
can result in large backgrounds for xenon-based dark matter experiments. The 
demonstrated background discrimination power of 1000 via pulse shape 
analysis is not sufficient to discriminate against electronic recoil events in- 
duced by /9-decay. However, it is likely that tritium will be reduced by a 
large factor during purification of the xenon. 

In Table[5]we tabulate the cosmogenic production rate of in various targets 
of dark matter detection experiments. These numbers can be used to guide the 
requirements for electron-recoil rejection and cryogenic purification to prevent 
background. 

It is worthwhile to mention that it has been demonstrated that liquid argon 
can achieve a pulse shape discrimination power of 1.3 x 10^ against electronic 
recoil events at energy of 15 keV [31]. Thus, at such a high energy threshold, 
the cosmogenic '^H will not be a problem for LAr-based dark matter detection 
and the dominant background source will be ^'^Ar contained in natural argon. 



5 Conclusions 

We have investigated the cosmogenic production of various isotopes in several 
target or source materials pertinent for dark matter and double-beta decay 
experiments. The tritium production in these materials due to cosmic-ray 
neutrons is substantial and steps must be taken to either reduce exposure 
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Table 4 

The production rate (per day per kg) of some long lived cosmogenic isotopes in 
natural xenon. Also shown is the half-life. Isotopes with very long half-lives or very 
small production rates are not listed because their contributions to the background 
are expected to be negligible. 



Cosmogenic Isotope 


Production Rate (/(kg day)) 


*l/2 




16.0 


12.33 y 


121ml 


11.7 


154 d 


123ml 


12.1 


119.7 d 


_LZ 1 [ILL fzi 


5.0 


109 d 




0.04 


3.3 y 




0.04 


2.7582 y 


119ml^^ 


0.02 


293.1 d 




0.004 


129.2 d 




3.2 


462.6 d 


113ml (27^ 


0.002 


14.1 y 



Table 5 

The cosmogenic production rate (per day per kg) of in various targets. 



Target 


Ar 


Xe 


Nal 


Csl 


Te02 


CaWOi 


Rate (/kg/day) 


44.4 


16.0 


31.1 


19.7 


43.7 


45.5 



of the target to cosmic rays, reduce the resultant '^H within the target after 
exposure, or develop an event-by-event analysis to remove decay events 
from the data stream. 
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